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We study the phase coherence and visibility of trapped atomic condensates on one-dimensional 
optical lattices, by means of quantum Monte-Carlo simulations. We obtain structures in the visibility 
similar to the kinks recently observed experimentally by Gerbier et al. [Phy. Rev. Lett. 95, 050404 
(2005); Phys. Rev. A 72, 053606 (2005)]. We examine these features in detail and offer a connection 
to the evolution of the density profiles as the depth of the lattice is increased. Our simulations reveal 
that as the interaction strength, U, is increased, the evolution of superfluid and Mott-insulating 
domains stall for finite intervals of U. The density profiles do not change with increasing U. We 
show here that in one dimension the visibility provides unequivocal signatures of the melting of 
Mott domains with densities larger than one. 
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The realization of trapped Bose-Einstein condensates 
in ultracold atoms on optical lattices has opened up the 
possibility of observing experimentally various quantum 
phases - e.g., superfluid (SF) and Mott-insulator (MI) - 
and the study of the nature of the transitions between 
them in a well-controlled manner. Indeed, the existence 
of SF and MI phases on optical lattices was established 
experimentally 0, Q, where it was demonstrated that 
by increasing the optical lattice depth the system passes 
from a SF phase to a predominantly MI one. Contrary to 
the unconfincd case, in traps there is in general a coexis- 
tence of SF and MI domains. Hence, the passage from SF 
to MI has to be understood as a cross-over rather than as 
a quantum phase transition 0, 0] , although a vestige of 
the latter remains in the guise of local quantum criticality 



The experimental systems can be modeled by the boson 
Hubbard model Q, described in one dimension (ID) by 



H = -t ( a ■ fli + i + af +1 flj j + fi rij 

i i 

+V T x i rii + U/lY.niirH-l), (1) 

i i 

where L is the number of sites and x% = ia is the coor- 
dinate of the ith site, and a is the lattice constant. The 
hopping parameter, t, sets the energy scale, m = a|<z, is 



Sij are bosonic creation 



the number operator, [aj,aj] 
and destruction operators. Vt is the curvature of the 
trap, while the repulsive contact interaction is given by 
U. The chemical potential, fi, controls the number of 
particles. The phase diagram of this model in the ab- 
sence of the confining trap has been extensively studied 
with the goal of elucidating the various quantum phases 
it exhibits H 0, and the transitions between them. 
The key experimental signature of these phases lies in 



the interference pattern observed after the release of the 
gas from the trap and subsequent free expansion - an SF 
(MI) produces a sharp (diffuse) interference pattern re- 
flecting the presence (loss) of phase coherence. Phase co- 
herence, especially in reduced dimensionality, continues 
to be of great interest both experimentally and theoret- 
ically. Particular attention has been focused recently on 
mechanisms which can destroy quasi-long range coher- 
ence in systems on optical lattices especially in ID 0. 
Our focus in this Letter is the role, in ID, of the passage 
from the SF to the MI phase in destroying phase coher- 
ence, which can be studied in matter wave interference. 

Whereas previous studies of SF-MI transition focused 
on the height Q and width Q of the central interfer- 
ence peak, an alternative scheme was proposed recently 
[llj where the reduction of phase coherence approaching 
the MI was characterized by the visibility of interference 
fringes, 
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Here S m&x and S m i n are the maximum and minimum val- 
ues of momentum distribution function, 
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It was observed that as the optical lattice depth [equiva- 
lent to the Hubbard U in Eq. QJ] is increased, the visibil- 
ity decreases until special values of U are reached where 
V displays "kinks" after which it decreases again. It was 
also shown that the values of U at which such kinks 
are observed are reproducible, and that they depend on 
the filling (number of atoms). A perturbative treatment 
of the homogeneous MI phase has shown that V de- 
creases as U~ 1 , improving on previous numerical results 
on small systems 
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Gerbicr et al. [ll| proposed that the kinks are linked 
to a redistribution of the density as the SF shells trans- 
form into MI regions with several atoms per site. In this 
Letter, we examine in detail the presence and properties 
of these visibility kinks with the help of quantum Monte 
Carlo (QMC) simulations of the boson Hubbard model 
using the Stochastic Series Expansion (SSE) method . 
We focus on ID optical lattices and show that while the 
kinks are indeed related to the redistribution in the den- 
sity associated with SF-MI transition, they are not solely 
produced by the transformation of SF shells into MI do- 
mains. Indeed, we find V reveals other subtle details of 
density redistributions with U. 

We start our study with the simplest case, i.e., a sys- 
tem in which the density in the middle of the trap never 
reaches n = 2, so that when the interaction strength is 
increased, only Mott domains with n = 1 appear. In Fig. 
H^a) we show the visibility and S max as functions of U/t. 
As in the experiments [ll|, V decreases with increasing 
U/t - reflecting the decrease of S max , and the increase 
of S m i n (not shown in the figure) - with an intermediate 
region over which it remains fairly constant. Two kinks 
can be observed both in V and S max . The first one (less 
evident) occurs around U/t — 6.1, and the second one 
around U/t = 7.0. 
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FIG. 1: (color online), (a) Visibility V and Smax as functions 
of the on-site interaction U/t. Initially, V decreases as U/t 
increases. After U/t ~ 6.3 its rate of reduction decreases 
due to the freezing of the density profiles (see text). The 
fast decrease after U/t « 7 is related to the formation of the 
central Mott core, (b) Density profiles at four different values 
of U/t, and in the Tonks- Girardeau (TG) regime. The profiles 
for U/t = 6.3 and 6.8 virtually coincide. The system under 
consideration has 40 bosons on a 80-site chain, and a trapping 
potential Vtcl 2 = O.Olt. Error bars on the data in this and all 
subsequent figures are smaller than the symbol sizes. 



Density profiles corresponding to four values oiU/t are 
depicted in Fig.^b). The density profile for U/t = 6.3 in 
Fig. ^b) shows that the first kink in Fig. ^a) (signaled 
by the first arrow) is related to the emergence of two 
MI plateaus at the sides (xi/a ~ ±(8 — 12)) of a central 
SF region. The second kink in Fig. ^a) is related to 
the formation of a full MI domain in the middle of the 
trap, which produces more evident structures in V and 
Smax- This occurs for U/t = 7.1 as shown in Fig. QJb), 
and signaled by the second arrow in Fig.^a). Plotting 
V and S m ax &s a function of U/t, allows us to present 
more precisely the position and shape of the kinks: In 
experiments, the control parameter is the ratio between 
the lattice depth and the recoil energy, which produces 

exponential changes in u/t 0,1m. 

One unexpected feature is the freezing of the density 
profiles before the full MI forms in the middle of the trap, 
which coincides with the plateau-like behavior of V and 
Smax between the two arrows. As U/t is increased be- 
tween 6.3 and 6.8 almost no changes occur in the density 
distribution, i.e., the bosons are no longer being pushed 
out of the central regions to the outlying zones even 
though U/t continues to increase. This behavior may 
seem surprising, as the central region of the system is SF, 
i.e., compressible, but can be explained by the presence 
of the emerging MI domains at the sides. The central SF 
region gets trapped between them, and the interaction 
U/t first has to increase a finite amount before particles 
can be transferred to the SF regions at the edge against 
the substantially larger trap energy there. 

This can be better understood by computing the total 
trapping (Ep) and interaction (Ep) energies as functions 
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FIG. 2: (color online), (a) Trapping (Et) and interaction 
(Ep) energies as functions of U/t. (b) Ratio 7 = \E P /Ek\ 
of potential to kinetic energy, and the chemical potential (fi) 
needed to maintain Nb = 40. The results are for the system 
of Fig. 



of U/t [Fig. In the interval U/t = 6.3-6.8 both 

quantities exhibit a plateau, which is also reflected in the 
chemical potential of the system [Fig. EJb)] . This occurs 
even though the total energy (not shown) increases con- 
tinuously, due to the continuous decrease in magnitude 
of the (negative) kinetic energy (Ek) of the bosons. One 
can then see that the formation of the full MI plateau 
is accompanied by a fast increase in the total trap en- 
ergy of the system by ~ At, the bandwidth in ID. On 
the other hand, the decrease of interaction energy pro- 
duced by the formation of the MI plateau is even larger 
~ 6t. Thus, in experiments, abrupt changes can occur in 
the density profiles even if the lattice depth is increased 
slowly. This can produce the escape of particles from the 
trap, heating, or other unexpected features. 

In Fig.Efb) we also show 7 = \Ep / Ek | , the ratio of po- 
tential to kinetic energy. This quantity is different from 
the one often used to characterize trapped bosons on lat- 
tices 7x, = U/t 14]. In contrast to jl, for the system in 
Figs. and 7 decreases with increasing U . This oc- 
curs because the density all over the trap becomes n < 1, 
and the double occupancy is strongly suppressed. (In the 
Tonks-Girardcau limit (TG), i.e., U — + 00, 7 = while 
7l = 00.) Like the visibility and the chemical poten- 
tial, 7 remains almost unchanged in the region where the 
density profiles are frozen. 

As the on-site interaction is further increased no more 
abrupt changes occur in the trap. The density profile 
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FIG. 3: (color online), (a) Visibility, V, as a function of 
the on-site interaction U/t, for Nb = 60, and Vra 2 = 0.06t, 
parameters which allow both n = 1 and n = 2 Mott regions to 
exist. For comparison, results for pure Mott insulating phases 
with n — 1 and n — 2 in open lattices without a trap are also 
given. In the inset the straight lines show the perturbative 
results of Ref. [TUl in ID (see text), (b) Integrated density 
over 20 lattice sites around the center of the trap. 
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FIG. 4: (color online). Density profiles corresponding to the 
points signaled by arrows in Figs. [3] and [K] The continuous 
line in (b) is the result in the TG regime. Nt, — 60 and 
V T a 2 = 0.06f. 

remains almost the same, as seen in Fig.^b), where we 
have also plotted the exact result in the TG limit. The 
visibility and S max , reduce continuously to Vtg — 0.39 
and S Z%_ = 1.3 (obtained using the approach presented 
111 Ref. |l5j). Notice that even when U — > 00 the visibility 
does not vanish, due to SF domains surrounding the MI. 

When the density at the center of the trap is higher, 
and exceeds two, the evolution of the visibility with the 
on-site repulsion exhibits an even richer structure. Re- 
sults for a system in that regime are presented in Fig. 
Ela). The visibility, up to U/t ~ 13, is very similar to 
Fig. d[a). Density profiles for three values of U in that 
interval are presented in Fig. 0][a). One can see that the 
emergence of MI regions with n — 1, and n = 2 surround- 
ing SF regions with 2 > n > 1, and n > 2, respectively, 
produces a plateau in V due to a freezing of the density 
profiles when increasing U. In Fig.^a), the formation of 
the n = 2 plateau abruptly reduces the visibility similar 
to the formation of the n = 1 plateau in Fig. da). 

However, the behavior above U/t = 13 has additional 
structures compared to Fig. da). In order to understand 
the origin of these visibility features, we have plotted 
in Fig. Efb) the integrated density over 20 lattice sites 
around the center of the trap N c = X^-io 71 *- ^ clear 
one to one mapping between the features in the visibility 
and plateau in N c is seen. The visibility kinks result not 
from the formation of new SF or MI regions, but rather 
from a redistribution of bosons between the MI states 
with n = 2, and n = 1. As seen in Figs.OJb), and0J such 
a redistribution occurs discontinuously in U. In addition, 
since the SF domains with 2 > n > 1 can increase their 
sizes during such a process, the visibility can increase [see 
for example the kinks around U/t = 14.6 and 21.5 in Fig. 
|3[a), and the corresponding density profiles in Fig. 0]. 

The above features are not restricted to the ID charac- 
ter of the system, and could be observed in higher dimen- 
sions. However, as U is increased even further (U > 25i), 
a purely ID effect sets in. As the MI plateau with n = 2 
melts, correlations start to develop between the two dis- 
connected SF domains with 2 > n > 1. This produces a 
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FIG. 5: (color online). Ratio 7 = \E P /Ek\ of potential to 
kinetic energy, and the total energy of the system (Es) vs U /t, 
for the system of Figs. EUandU N b = 60 and V T a 2 = 0.06t. 

large increase of the visibility, as seen in Fig. EH a) . [The 
corresponding density profiles are shown in Fig.^Jb).] In 
ID this increase in the visibility provides an unambiguous 
signature of the presence, and melting, of the n = 2 (or 
larger) MI domain. This can be useful for understanding 
the dynamics of strongly correlated bosons in ID |l6l| . 

For very large values of U, beyond the ones in Fig. 
OH a), the 2 > n > 1 SF domain will eventually disappear, 
as occurs in Fig, ffla). producing a further reduction in 
the visibility. In the TG regime, we obtain (for these pa- 
rameters) Vtg — 0.02. The corresponding density profile 
can be also seen in Fig.HJb). 

We have also plotted in Fig. E[a), the values obtained 
for the visibility in homogeneous systems with 60 bosons 
and densities n = 1 and n — 2. (We have used open 
boundary conditions as they are the closer to the trapped 
case) . These results in homogeneous systems are very dif- 
ferent from the ones in the trapped case. Due to the ex- 
istence of SF domains, the visibility in the trap is always 
larger than that in the homogeneous case. In the region 
of interest, where the MI plateau emerges, and melts, no 
extrapolation is possible from the uniform case. Only for 
very large values of U, after a MI domain appears in the 
center of the trap, can one can expect the uniform and 
trapped systems to behave similarly. In the inset we have 
compared the results for the homogeneous systems, with 
those obtained in Ref. HH [Vld = 4(n + l)t/U]. For the 
largest values of U one can see that the t/U power law 
starts to develop, but its prefactor is still different from 
4(n +1), so that very large values of U are needed for a 
good agreement in ID. 

We conclude by showing in Fig. [3] the behavior of 7 in 
the system of Figs. [3] and 0] In this case, since the density 
at the middle of the trap is larger than one, i.e., there is 
significant double occupancy in this region, 7 (and Ep) 
increases with U/t. It also exhibits the same jumps pro- 
duced by the redistribution of particles in Fig. 0] As in 
the system in Figs. ^ an( l HI this occurs even when the 
total energy of the system (Es) increases continuously 
with U, as can be also seen in Fig. 

In this Letter we have explored the evolution of the 



visibility of trapped atomic gases in one dimensional op- 
tical lattices using Quantum Monte Carlo simulations. 
We have shown that the visibility behaves very similar 
to that observed experimentally. In particular, it has 
kinks associated with redistribution of density amongst 
Mott insulating and superfluid regions within the trap. 
In addition, we have also exhibited several other novel 
features of the visibility evolution in ID, like a large in- 
crease due to the melting, with increasing U/t, of n > 1 
MI plateaus. We have demonstrated that the evolution 
of the density distribution with interaction strength ex- 
hibits pauses. That is, at certain values of U the density 
distribution, and other obscrvables, do not change even 
when the interaction strength increases over a range as 
large as t/2. We have shown that the emergence of this 
static behavior is associated with the formation of Mott 
insulating plateaus away from the trap center. These 
plateaus block the transfer of bosons to the outer parts 
of the system, and hence cause the evolution to stall. 
While many quantities in trapped Bose systems are well 
described by the local density approximation, it is not 
clear that approach will capture the above behavior, in- 
cluding the kinks in the visibility. This is because these 
effects arc intrinsically tied to the competition of the trap 
versus kinetic and interaction energies in systems where 
the SF and MI domains are of finite width, as in the ones 
explored in the recent experiments. 
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